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This  report  results  from  a  contract  tasking  University  of  BARI  as  follows:  The  kinetics  of  nozzle  expansion  flow  is  a  topic  of  large  interest 
due  to  its  linking  with  technological  applications.  Usually  one  can  treat  the  problem  by  using  two  approaches.  The  first  one  is  to  use 
complicated  2D  fluid  dynamics  codes  with  poor  kinetics,  the  second  one  is  to  use  ID  nozzle  codes  emphasizing  the  role  of  chemical 
kinetics  in  affecting  the  whole  fluid-chemical  problem.  The  second  approach  has  been  used  by  our  group  to  characterize  simple  gases  such 
as  F12,  N2,  02  and  their  mixtures.  The  peculiarity  of  our  approach  is  to  describe  each  vibrational  level  of  the  molecule  as  a  new  species 
thus  avoiding  the  concept  of  vibrational  temperature.  This  approach  could  be  in  principle  extended  to  polyatomic  molecules  even  though  in 
this  case  the  state  tot  state  kinetics  becomes  a  very  difficult  problem.  These  concepts  will  be  developed  in  the  present  proposal  which  deals 
with  the  characterization  of  the  kinetics  of  NF13  through  nozzles  of  different  geometries.  Our  proposal  consists  of  different  steps  including: 

1 )  development  of  a  realistic  kinetics  for  ammonia  decomposition,  2)  development  of  a  realistic  macroscopic  model  for  vibrational 
deactivation  of  polyatomic  molecules,  3)  evaluation  of  the  dependence  of  rates  on  vibrational  temperature,  4)  insertion  of  points  1-3  in  the 
nozzle  equations  and  5)  examination  of  results  for  different  nozzle  geometries. 
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Introduction 

Recently,  solar  thermal  propulsion  (STP)  engines  have  been  investigated  as  a 
promising  technology  to  improve  the  capabilities  of  small  satellites.  In  the  Surrey 
Space  Centre  [1,2],  a  microscale  engine  has  been  designed  to  use  ammonia  as 
propellent.  The  main  advantages  of  ammonia  propellent  are  its  sufficiently  low  molar 
mass  and  its  chemical  inertia  at  room  conditions  (not  as  dangerous  as  hydrazine). 

To  get  better  performances  the  ammonia  gas,  in  the  heat  exchanger,  should  be 
dissociated  in  N2  and  H2  species,  doubling  the  particle  density.  Unfortunately  the 
dissociation  kinetics,  at  temperatures  of  interest  in  the  engine,  is  so  slow  that  the 
complete  thermal  dissociation  of  ammonia  molecules  could  not  occur. 

In  this  project  we  have  developed  a  kinetic  model  for  an  ammonia  mixture  flowing 
into  the  nozzle  investigated  in  previous  studies  [1,2].  It  is  worth  noting  that  we  have 
introduced  a  state-to-state  kinetics  to  verify  the  presence  of  non-equilibrium 
vibrational  distributions  during  the  nozzle  expansion. 

Nozzle  Geometry 

The  nozzle  used  in  the  thruster  has  a  conic  geometry  of  the  following  dimensions  (in 
mm): 


throat  radius  0.35 

exit  radius  2.0 

exit-throat  distance  6.2 

inlet  radius  1.0 

inlet-throat  distance  1 .2 


These  data  have  been  provided  by  the  Surrey  Laboratories.  The  converging  nozzle  is 
quite  short  so  we  choose  to  include  part  of  the  reservoir  in  the  calculation  through  a 
connecting  cone  of  the  following  dimensions  (in  mm) 

reservoir  radius  1 1 .25 

connecting  cone  length  2.0 

The  nozzle  is  presented  on  fig  1 . 

The  fluid  dynamic  model  of  the  nozzle  flow  is  based  on  the  almost  one-dimensional 
Euler  equations  for  a  variable  section  duct.  The  Euler  equations  are  solved,  at  the 
stationary  state,  together  with  the  perfect  gas  equation  of  state  and  the  mass  continuity 
for  each  species  including  the  reaction  contributions.  The  numerical  scheme  adopted 
is  described  in  [3], 

Chemical  Model 

To  build  the  chemical  kinetic  models,  we  need  to  select  both  the  species  that  must  be 
taken  into  account  and  their  mutual  reactions.  In  the  laboratory  experiments,  the 
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pressure  ranges  from  1  to  4  atm  and  the  temperature  is  less  than  1500  K.  The  pressure 
values  are  measured  directly,  but  the  gas  temperature  is  estimated  only  from  some 
measurements  performed  outside  the  reservoir.  This  results  in  a  large  incertitude  on 
the  actual  conditions  operating  inside  the  reservoir.  In  these  context  the  ammonia 
decomposition  kinetics  is  quite  slow.  Surface  dissociation,  neglected  in  this 
preliminary  approach  should,  probably,  be  inserted  in  the  model. 


Fig.  1:  Nozzle  profile. 


x(m) 


1.  Species 

We  assume  that  the  ammonia  mixture  consists  of  atoms  and  molecules  formed  in  the 
NH3  decomposition.  We  ignore  all  the  ions  because  the  temperatures  involved  are 
quite  low.  We  introduce  the  vibrational  excitation  only  for  N2  and  H2  molecules, 
which  have  stable  vibrational  excited  levels.  For  all  the  other  molecules,  we  diregard 
the  internal  levels  as  their  number  is  so  high  that  the  computational  load  becomes 
unaffordable.  Multitemperature  model  could  be  a  valid  approach  assuring  an 
acceptable  computational  effort  nevertheless,  in  this  work,  we  consider  these  internal 
states  in  equilibrium  with  the  translational  degree  of  freedom  or  as  following  a  frozen 
kinetics.  According  to  the  internal  state  model  selected  we  assign  the  Cp  and  y  (y  = 
Cp/Cy)  coefficients  to  each  molecule.  For  each  vibrational  mode  or  rotational  axis,  in 
equilibrium  with  the  translational  degree  of  freedom,  we  add  the  quantity  aR  into  the 
expressions  of  Cp  and  Cv,  where  R  represents  the  universal  gas  constant  and  the  a 
parameter  is  set  to  1  for  a  vibrational  mode  and  to  !4  for  a  rotational  axis.  This 
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approximation  is  valid  for  the  rigid  rotor  and  for  the  infinite  harmonic  oscillator,  since 
in  general  the  Cp  coefficient  depends  on  the  temperature.  For  atomic  species  no 
internal  degree  of  freedom  is  considered.  The  contribution  of  translational  degrees  to 
the  molar  specific  heat  is  1.5R.  Table  1  reports  the  species  accounted  and  their 
relative  properties. 


Tab.  1:  species  inserted  and  their  physical  properties.  Rot  is  the  number  of 
independent  rotational  axis,  Vib  the  number  of  vibrational  modes  and  Ef  the  enegy  of 
formation. 


Species 

Mass  (AMU) 

Rot 

Vib 

Ef(eV) 

nh3 

17 

3 

6 

-0.403242 

N2H3 

31 

3 

9 

4.30994 

nh2 

16 

3 

3 

2.00285 

NH 

15 

2 

1 

3.90225 

n2 

28 

2 

1 

0.00000 

h2 

2 

2 

1 

0.00000 

N 

14 

0 

0 

4.87950 

H 

1 

0 

0 

2.23910 

For  H2  and  N2  we  estimate  the  energy  of  each  vibrational  level  (ev)  through  the  semi- 
empirical  formula,  function  of  the  vibrational  quantum  number  and  of  some 
spectroscopic  constants  expressed  in  cm"1: 


s 


v 


(DgXg 


+  (ney  g 


+  (DgZg 


The  H2  spectroscopic  constants  have  been  obtained  through  the  fits  of  the  vibrational 
level  energies  extracted  from  the  BKMP  potential  energy  surface  [4]  while  for  N2  we 
used  data  from  ref.  5.  Table  2  summarizes  those  values  for  the  two  species. 

Tab.  2:  Spectroscopic  constants  (in  cm"1)  used  in  the  calculation  of  the  vibrational 
level  energies. _ _ _ 


Parameter 

n2 

h2 

0)g 

2358.57 

4464.7 

(OgXg 

14.324 

146.93 

weye 

-0.00226 

4.7108 

(DgZg 

0.0 

-0.22712 

Ediss 

9.7639 

4.4772 

V  max 

48 

15 

2.  Processes 

We  implement  in  our  model  a  macroscopic  kinetics  for  NH3  depletion  introducing  the 
following  processes: 


NH3+X  NH2+H+X 

nh3  +  X  nh  +  h2  +  X 
nh3  +  h  nh2  +  h2 


(pi) 

(p2) 

(p3) 
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NH  +  H  +  X  ^  NH2  +  X 

(p4) 

H2  +  N  +  X  NH2  +  X 

(p5) 

NH2  +  H  NH  +  H2 

(p6) 

NH2  +  NH2  NH3  +  NH 

(p7) 

N  +  H  +  X  NH  +  X 

(p8) 

NH  +  H  H2  +  N 

(p9) 

NH  +  N  N2  +  H 

(plO) 

NH  +  NH  H2  +  N2 

(pH) 

2H  +  N2  H2  +  N2 

(pl  2) 

2H  +  H2  H2  +  H2 

(pl  3) 

2N  +  N2  °  N2  +  N2 

(pl  4) 

2N  +  N  N2  +  N 

(Pl5) 

NH3  +  NH2  -*•  N2H3  +  h2 

(p!6) 

N2H3  +  X  ->  NH2  +  NH  +  X 

(Pl7) 

N2H3  +  H  -*  2NH2 

(P 1 8) 

N2H3  +  H  ->  NH3  +  NH 

(p!9) 

where  X  is  a  generic  component.  The  rate  coefficients  for  all  these  processes  are 
expressed  in  the  Arrhenius  form: 

K  =  K0Tae-Ea/T 

For  the  direct  reaction,  the  parameters  Ko,  a  and  Ea  are  reported  in  table  3.  The  rate 
coefficients  of  these  processes  have  been  taken  from  ref.  [5-10].  In  this  model,  we 
neglect  more  complex  molecules,  such  as  N2H4,  that  do  not  seem  to  be  significant. 
Rates  for  the  reverse  processes  have  been  obtained  through  the  detailed  balance 
principle 


Keq 

where  Keq  is  the  equilibrium  constant.  Equilibrium  constants  have  been  calculated, 
from  partition  functions,  following  the  classical  statistical  thermodynamic  theory  [11] 
subsequently  they  have  been  fitted  by  the  general  equation 

T 

ln(Keq)  =  K„+Kpf-^j 

We  consider  only  the  equilibrium  constants  for  the  generation  of  a  new  species  from 
standard  ones  like  N2  and  H2  in  nitrogen-hydrogen  compounds  formation.  Therefore 
all  the  formation  reactions  can  be  synthesized  as 

fN2  +  |H2-NnHm 

Table  4  presents  the  parameter  values  utilized  to  calculate  Keq,  for  the  species  listed  in 
table  1. 
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Tab.  3:  Arrhenius  coefficients  for  the  reactions  inserted  in  the  model.  Ko  has  the 
dimension  of  cnr  (n_1Vs,  where  n  is  the  number  of  particles  involved  in  the  reaction,  Ea 
is  expressed  in  K  and  a  is  a  dimensionless  parameter. 


process 

Ko 

a 

Ea 

Ref. 

Pi 

1.53e-8 

0 

42400 

[6] 

p2 

1.00e-9 

0 

47000 

T4] 

p3 

3.20e-13 

0.67 

1720 

T5] 

p4 

1.00e-32 

0 

0 

[7] 

p5 

1.00e-34 

0 

0 

[7] 

p6 

2.30e-13 

0.67 

2160 

T51 

P7 

6.64e-12 

0 

2800 

T61 

P8 

1.00e-33 

0 

0 

m 

p9 

1.70e-12 

0.68 

957 

T5] 

plO 

1.77e-l  1 

0 

0 

[5] 

pH 

6.60e-13 

0.55 

957 

T5] 

pl2 

1.34e-31 

-0.60 

0 

T5] 

pl3 

2.68e-31 

-0.60 

0 

[5] 

pl4 

7.44e-32 

-0.50 

0 

[5] 

p  1 5 

3.31e-27 

-1.50 

0 

T51 

pl6 

1.33e-12 

0.5 

10850 

T4] 

Pl7 

1.70e-8 

0. 

21000 

[4] 

pl8 

2.60e-12 

0 

0 

T4] 

pl9 

1.70e-13 

0 

0 

[4] 

Tab  4:  Coefficients  used  to  calculate  the  formation  equilibrium  constants  when  the 
pressures  are  expressed  in  Pa. _ _ _ _ 


Species 

Koo 

Kp 

9p 

Ke 

9e 

nh3 

51.8293 

-13.1631 

0.960208 

-7.97367 

628.095 

nh2 

19.6666 

63.4542 

1.01863 

13.5452 

311.287 

NH 

-5.14978 

85.2911 

0.998065 

-2.21874 

423.673 

N 

28.3772 

-114.298 

0.997620 

-2.14877 

4371.66 

H 

26.1506 

-52.7759 

0.996362 

-2.71911 

2157.34 

From  these  formation  equilibrium  constants  we  can  calculate  the  equilibrium 
constants  relative  to  all  the  reactions  involving  the  species  in  the  table. 

The  kinetic  model  has  been  improved  including  a  state-to-state  kinetics  for  the 
diatomic  species  N2  an  H2,  substituting  processes  (P12-P15)  with  the  state-selective 
dissociation  and  adding  vibrational  relaxation  processes  [12-14].  All  the  other 
processes  involve  only  the  ground  state  of  such  molecules. 

Results 

The  kinetics  of  ammonia  is  not  well  known,  even  if  this  subject  has  been  widely 
studied.  The  main  problem  is  the  kinetics  slowness  that  prevents  from  estimating  the 
kinetic  constant  precisely.  Experimental  works  [1,2]  on  ammonia  thrusters  assume  a 
complete  decomposition  of  ammonia  in  the  reservoir,  but  the  kinetic  scheme  proposed 
does  not  permit  to  reach  the  equilibrium  condition  in  an  acceptable  time  for  T<2000 
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K,  even  if  the  asymptotic  value  agrees  to  the  equilibrium  composition.  For  this  reason 
we  compare  the  results  obtained  starting  from  pure  ammonia  and  from  equilibrium 
concentrations. 


V) 

a 

© 

'-G 

© 

a 

j 

~o 


P  =4  Bar  ;  T=1500K 

o  o 


inlet  throat 


x(m) 

Fig.  2:  Molar  fractions  of  mixture  species  along  the  nozzle.  We  assume  as  reservoir 
conditions  Po=4  Bar  and  To=1500  K  and  start  from  an  equilibrium  inlet  composition. 

We  investigate  four  different  pressure  values  in  the  range  1-4  Bar.  In  figure  2  we 
report  the  concentration  profile  along  the  nozzle  obtained  for  reservoir  conditions  of 
P=4  Bar  and  T=1500  K,  and  introducing  the  vibrational  kinetics  of  N2  and  Ffi. 

Tab.  5:  Mass  flow  rate  and  thrust  of  the  nozzle  as  a  function  of  the  inlet  pressure.  The 
inlet  temperature  is  To=1500  K  except  when  marked  with  the  “*”  symbol  assigned  to 
T0=2500  K.  At  inlet  we  suppose  an  equilibrium  composition. 


Po(Bar) 

mass  flow  rate  (g/s) 

thrust  (mN) 

1 

2.1708e-2 

64.337 

2 

4.3404e-2 

128.68 

3 

6.5077e-2 

192.97 

4 

8.6784e-2 

257.38 

4* 

6.5954e-2 

256.74 

*  T=2500  K 
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5  10 

vibrational  quantum  number 

Fig.  3:  Hydrogen  vibrational  distributions  at  different  nozzle  positions  in  the  same 
conditions  as  in  fig.  2. 


Tab.  6:  Mass  flow  rate  and  thrust  of  the  nozzle  determined  experimentally  in  different 
test  cases. 


mass  flow  rate  (g/s) 

thrust  (mN) 

.140 

219 

.147 

223 

.117 

183 

.121 

190 

Tab  7:  Mass  flow  rate  and  thrust  as  a  function  of  the  inlet  dissociation  degree  (a)  of 
the  ammonia.  Inlet  pressure  and  temperature  are  Po=4  Bar  and  To=1500  K. 

a 

mass  flow  rate  (g/s) 

thrust  (mN) 

0.0 

1.11 89e- 1 

284.13 

0.2 

1.0437e-l 

280.35 

0.4 

9.7946e-2 

275.67 

0.8 

8.9483e-2 

264.37 

1.0 

8.5450e-2 

253.33 
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N  VDF 


—  x(m)=-0.0058 
*■  -  -x(m)==-0.0012 
•-  -x(m)==0 

- x(m)==0.0001 

. x(m)==0.001 

—  x(m)==0.006 


vibrational  quantum  number 

Fig.  4:  Nitrogen  vibrational  distributions  at  different  nozzle  positions  in  the  same 
conditions  as  in  fig.  2. 
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T  =1500 K 

o 


inlet  throat 


x(m) 

Fig.  5:  Pressure  profile  along  the  nozzle  as  a  function  of  the  inlet  pressure  for 
To=1500  K  and  starting  from  an  inlet  equilibrium  composition. 

It  is  interesting  to  observe  that  the  flow  is  practically  frozen  along  the  whole  nozzle 
except  for  some  minor  species  such  as  H  and  NH2  that  completely  disappear  at  the 
nozzle  exit.  The  atomic  nitrogen  molar  fraction  is  lower  than  10"10. 

Figures  3  and  4  describe  the  vibrational  distribution  of  hydrogen  and  nitrogen  for  the 
conditions  of  figure  2.  We  can  observe  that  the  nitrogen  distributions  are  frozen.  On 
the  contrary,  while  Ffi  distributions,  for  v<3,  are  almost  frozen,  their  tails  are  strongly 
depleted,  resulting  in  non-Boltzmann  distributions.  The  different  behaviour  of  the  two 
species  is  due  to  the  different  percentage  of  nitrogen  and  hydrogen  atoms;  hydrogen 
atom  concentration  is  high  enough  so  that  the  VT  collisions  involving  atoms  are  fast 
enough  to  cool  the  distribution  tails.  We  must  point  out  that  we  neglect  the  vibrational 
depletion  of  N2(v)  by  atomic  and  molecular  hydrogen  collisions.  We  observe  that  the 
Mach  number  is  practically  independent  on  the  initial  temperature  and  pressure, 
obtaining  at  the  nozzle  exit  a  value  close  to  Mach  5.4.  We  report  the  pressure  along 
the  nozzle  in  figure  5.  It  is  quite  interesting  to  note  that  the  exit  pressure  is 
proportional  to  the  inlet  pressure.  These  two  results  occur  as  the  chemistry  weakly 
affects  the  macroscopic  quantities  during  the  expansion. 
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Table  5  displays  the  mass  flow  rate  and  the  thrust  for  different  inlet  pressures  at 
To=1500  K  or  To=2500  K.  The  agreement  with  the  experimental  results  reported  in 
table  6  is  satisfactory,  as  the  nozzle  reservoir  conditions  (mainly  the  temperature)  are 
not  well  known  experimentally. 

There  are  some  incertitude  in  the  chemistry  model,  because  the  dissociation  of 
ammonia  is  very  slow  for  T<2000  K.  In  such  conditions  the  kinetic  models  published 
in  the  literature  cannot  provide  the  complete  dissociation  of  ammonia  in  an  acceptable 
time  period.  On  account  of  this  we  also  try  as  initial  condition  a  pure  ammonia  gas. 
Dissociation  inside  the  reservoir  and  the  nozzle  is  very  poor,  so  that  again  in  this  case 
the  flow  is  frozen.  As  a  consequence,  the  flow  properties  depend  only  on  the  Cp/Cv 
ratio  and  on  the  mean  molar  mass,  then  the  exit  properties  are  completely  different 
form  the  equilibrium  case.  As  an  example  in  figure  6  we  compare  the  Mach  numbers 
computed  starting  from  an  equilibrium  composition  (Eq.)  and  from  pure  ammonia 
(NoEq.).  It  is  not  possible  to  really  know  what  composition  prevails  inside  the 
reservoir  for  T<2000K,  moreover  some  impurites  or  wall  catalitic  effects  can 
accelerate  the  kinetics,  bringing  the  system  to  equilibrium  faster  than  predicted  by  our 
kinetic  model. 
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Fig.  6:  Mach  number  profile  compared  for  inlet  equilibrium  (Eq)  and  pure  ammonia 
(noEq)  case  (pressure  and  temperature  conditions  as  in  fig  2). 

We  try  the  same  comparison  at  To=2500  K.  In  this  case,  the  mixture  reaches  the 
equilibrium  composition  in  the  reservoir.  Nevertheless  the  Mach  number  profiles  in 
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the  equilibrium  and  non  equilibrium  case  are  again  different  (see  fig.  7).  This  is  the 
demonstration  that  the  history  of  the  gas  in  the  reservoir  is  important. 
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Fig.  7:  Same  as  in  fig  6  with  reservoir  temperature  To=2500  K. 

Effects  of  initial  conditions  propagate  also  on  the  temperature  profiles.  In  figure  8  we 
show  the  normalized  temperature  profile  calculated  starting  from  an  equilibrium 
condition  (Eq.)  and  from  an  undissociated  ammonia  (NoEq.)  for  a  reservoir 
temperature  of  1500  K  and  a  pressure  of  4  Bar.  We  observe  that  the  exit  temperature 
is  much  higher  in  the  case  of  the  equilibrium  condition.  This  behaviour  could  be  due 
to  the  lower  mean  molar  mass  associated  to  the  mixture  in  equilibrium  composition 
inside  the  reservoir. 

Different  behaviours  occur  when  the  reservoir  temperature  reaches  2500  K  (see  fig. 
9).  In  fact,  at  this  temperature,  the  reactions  in  the  inlet  are  faster  than  in  the  flow, 
therefore  one  can  observe  an  increase  of  the  temperature  due  to  the  mixture 
dissociation.  Even  if  the  ammonia  energy  is  lower  than  that  of  nitrogen  and  hydrogen 
molecules,  the  compression  due  to  the  double  particle  number  density  produces  an 
increase  of  the  temperature.  At  the  nozzle  exit,  the  differences  are  lower  than  in  the 
case  of  an  inlet  temperature  of  1500  K. 

We  find  similar  behaviours  for  the  relaxation  of  the  vibrational  temperature, 
calculated  as  the  temperature  of  the  first  excited  level  (figure  10).  Starting  from  non 
equilibrium  conditions  we  can  observe  an  increase  of  the  vibrational  temperatures 
following  the  gas  temperature.  At  the  nozzle  exit  the  vibrational  temperature  freezes 
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at  a  temperature  higher  than  2300  K  just  after  the  throat.  N2  vibrational  temperature  is 
practically  frozen. 

As  a  consequence  of  the  previous  results,  we  can  state  that  there  are  some  incertitudes 
about  the  inlet  composition.  At  this  purpose  we  have  reported  in  table  7  the  calculated 
mass  flow  rate  and  thrust  for  different  dissociation  degree  “a”  of  ammonia  for  the 


a  =  1  - 


in 

II 

X 

z 

h2+ 

[NA 

[NH3]  +  ,5 

n2] 

that  ranges  from  0  (pure  ammonia)  to  1  (N2  and  H2  mixture).  We  note  that  both 
quantities  decrease  when  the  dissociation  degree  increases,  but  the  thrust  changes 
more  slowly  than  the  mass  flow  rate. 
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Fig.  8:  Reduced  temperature  profde  compared  for  inlet  equilibrium  (Eq)  or  pure 
ammonia  (noEq)  case,  for  a  reservoir  temperature  To  =  1500K  (pressure  and 
temperature  conditions  as  in  fig  2). 

Conclusions. 
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In  this  study  we  investigate  the  role  of  chemical  kinetics  in  modeling  ammonia  solar 
thermal  engine.  Many  problems  still  need  to  be  solved.  In  particular,  the  kinetic 
model  has  to  be  improved,  insertion  of  some  impurities  or  wall  catalysis  could  speed 
up  the  dissociation  kinetics. 

A  further  improvement  of  the  model  could  be  the  inclusion  of  non  equilibrium  in 
rotational  and  vibrational  levels  of  species  other  than  N2  and  H2,  even  if  it  seems  that, 
in  the  conditions  the  thruster  runs,  the  vibrational  kinetics  is  not  so  much  effective. 


T  =2500  K 

0 


x(m) 

Fig.  9:  Reduced  temperature  profde  compared  for  inlet  equilibrium  (Eq)  or  pure 
ammonia  (noEq)  case,  for  a  reservoir  temperature  To  =  2500K  (pressure  and 
temperature  conditions  as  in  fig  2). 
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T  =2500  K  P  =4  Bar 


x(m) 

Fig.  10:  Vibrational  temperature  (Tv)  profile  of  N2  and  H2  in  the  same 

conditions  as  in  figure  9. 
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